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Abstract. We investigate the interaction of thermonuclear flames in Type la supernova explosions with vortical flows by means 
of numerical simulations. In our study, we focus on small scales, where the flame propagation is no longer dominated by the 
turbulent cascade originating from large-scale effects. Here, the flame propagation proceeds in the cellular burning regime, 
resulting from a balance between the Landau-Darrieus instability and its nonlinear stabilization l Ropke et al., 2003b g). The 
interaction of a cellularly stabilized flame front with a vortical fuel flow is explored applying a variety of fuel densities and 
strengths of the velocity fluctuations. We find that the vortical flow can break up the cellular flame structure if it is sufficiently 
strong. In this case the flame structure adapts to the imprinted flow field. The transition from the cellularly stabilized front to the 
flame structure dominated by vortices of the flow proceeds in a smooth way. The implications of the results of our simulations 
for Type la Supernova explosion models are discussed. 
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1. Introduction 

In the canonical astrophysical model Type la supernovae (SNe 
la in the following) are associated with thermonuclear ex - 
plosions of white dwarf (WD) stars faovle & Fowlerl [i960). 
These WDs are composed of carbon and oxygen. Due to the 
high temperature sensitivity of the activation energy of carbon 
and oxygen burning, the reaction is confined to a very narrow 
region and propagates as a combustion wave. This defines a 
thermonuclear flame. The theoretical understanding of flame 
propagation in WDs provides the key to SN la explosion mod- 
els with the effective flame velocity being the most important 
parameter. What are the phenomena that mediate the combus- 
tion wave? The answer to that question is given by the hydro- 
dynamics of combustion. Conservation laws across the flame 
admit two fundamental modes of flame propagation: the det- 
onation, in which the flame propagates due to shock waves, 
and the deflagration, which is based on microphysical trans- 
port phenomena such as heat conduction and species diffusion. 

Flame propagation in the detonation mode is certainly eas- 
iest to model. Here the flame propagation velocity with respect 
to the unburnt material is larger than the corresponding sound 
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speed. Therefore the WD star has no time to expand prior to 
the incineration in this detonation model of SN la explosions 
and all the material is consumed at high densiti es resulting i n 
ashes consisting solely of iron-peak elements (lArnettLll969h . 
However, this is in conflict with observational spectra, which 
show strong indication of intermediate mass elements. 

Consequently, the star has to expand considerably before 
it is processed by the thermonuclear flame. This is possible 
if flame propagation starts out in the deflagration mode, since 
here the microphysical transport processes lead to a subsonic 
flame velocity. In the present study, we refer to that deflagra- 
tion model of SNe la. Moreover, underlying to our model is 
the so-called single-degenerate scenario, where the progeni- 
tor is a binary system with the WD accreting matter from the 
non-degenerate companion until it reaches the Chandrasekhar 
mass. At this point, thermonuclear reaction ignites near the 
center of the WD and propagates outward as a deflagration 
flam e. For a review of SN la explosion models we refer to 
iHillebrandt & NiemevelfcOOOl) . 

In the framework of the deflagration model, recent large- 
scale SN la explosion models have been quite successful 
jReinecke etaill2002tlGamezo et alll2003l) . Nevertheless, the 
determination of the flame propagation speed is much harder 
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in this case. A planar deflagration flame would be far too 
slow to cause powerful SN la explosions, but these mod- 
els predict an acceleration of the deflagration flame due to 
interaction with turbulence. The resulting energy release is 
suffi cient to unbind the WD st ar. Following the reasoning 
by iNiemever & Wooslevl dl997h . the turbulent motions are 
evoked by large-scale instabilities — such as the Rayleigh- 
Taylor (buoyancy) instability and the Kelvin-Helmholtz (shear) 
instability — a nd decay to smalle r scales forming a turbulent 
eddy cascade ( Richardsonl 1 1 922l) . Eddies of that cascade wrin- 
kle the flame and increase its surface. This corresponds to an 
acceleration of the effective flame propagation speed. However, 
the vast range of spatial scales relevant to the problem of SN la 
explosions (from the radius of the star down to the width of a 
typical flame it covers about 1 1 orders of magnitude) restricts 
these models to the simulation of effects on the largest scales 
only. This makes it impossible to directly determine the turbu- 
lent flame propagation velocity in the simulations and therefore 
assumptions on the physics on unresolved small scales must be 
made. The fundamental idea here is that flame propagation is 
dominated by the turbulent cascade that is driven solely from 
large-scale instabilities and that no other effects contribute to 
the generation of turbulence. 

However, due to the scaling law of the turbulent cascade, 
which predicts the turbulent velocity fluctuations to monoton- 
ically decrease with smaller scale, there exists a certain length 
scale below which the flame burns faster through the corre- 
sponding eddies than these can deform it. Fo r chemi cal flames, 
this scale is identified with the Gibson scale dPetersl ll986) and 
we will use this term in the context of thermonuclear combus- 
tion, too. Below the Gibson scale the flame propagates in a 
"frozen turbulence field" and is no longer affected by the eddy 
cascade. Large scale models assume stable flame propagation 
here. Although some theoretical concepts support this idea, no 
hydrodynamical simulations provided convincing proofs yet. 

In order to fill that gap in SN la models we performed 
numerical studies of flame propagation around the Gibson 
scale. The basic physical effects that determine flame propaga- 
tion here are the hydrodynamic al Landau-Darrieus instability 
(Darrieus, 1938; Landau, 1944) and the counteracting nonlin- 
ear stabilization of the flame dZerdovichl 11964 . The balance 
of the two effects gives rise to a stable cellular flame structure 
characterizing the cellular burning regime. The goals of our in- 
vestigations were to determine whether or not the cellular burn- 
ing regime exists for thermonuclear flames in WD matter and 
to test how robust the anticipated cellular stabilization is un- 
der varying conditions. On the one hand, these studies examine 
a fundamental assumption of present SN la explosion models 
and on the other hand they are intended to explore the possi- 
bility of new physical effects that may have to be included into 
these. In particular, the question of the robustness of the cel- 
lular stabilization of the flame is closely connected to the so- 
called delayed detonation models of SN la explosions. In these 
models, the flame starts out as a deflagration wave, but later 
turns into a supersonic detonation. Empirical one-dimensional 
SN la models support this idea, since those lead to quite real- 
istic spectra if a deflagration-to-detonati on trans ition (D DT) is 
imposed artificially at low fuel densities faoflich & Khokhlovt 



Il996l llwamoto et ail Il999h . Although DDTs are frequently 
observed in terrestrial combustion, a mechanism that could ac- 
count for such a transition under conditions of SN la explo- 
sions could not be identified so far jNiemeverlfl9 99). However, 
one possibility is a breakdown of the flame stabilization and 
a subsequent generation of turbulence on small scales by the 
flame itself, that would result in an additional acceleration of 
the flame. This active turbulent combustion (ATC) mechanism 
has been proposed by Niemever & Wooslev ( 1997) on t he ba- 
sis of a work by Kerstein ( 1996). Niemever & Wooslev dl997l) 
suggest that the cellular flame stabilization could break down at 
the Gibson scale leading to ATC. One motivation for our study 
is to examine this idea. 

The present paper is the third in a series of publications 
on the investigation of the cellular burning regime in Type la 
supernova explosions. After the introduction of the numerical 
methods to study this r egime together with providing some ex- 
amples of application ( Rop ke et all l2003bl) . the implementa- 
tion was used to stu dy flame propagation into quiescent fuel 
(Ropk e et all l2003al) . The results obtained there (to be sum- 
marized in the next section) raised the question of the interac- 
tion of a cellular flame with turbulent flows. This topic is also 
motivated from the supernova model. Around the Gibson scale 
one can still expect turbulent velocity fluctuations of the or- 
der of the laminar burning speed of the flame originating from 
the turbulent cascade. Additionally, pre-ignition convection in 
the WD may result in relic turb ulent motions of considerable 
strength. Hoflich & Stein ( 2002) claim surprisingly high values 
of convective velocities. 

A brief description of the physical effects underlying the 
cellular burning regime and of the numerical methods applied 
to model it will be provided in the next section. In Sect.[3]the 
results of our numerical simulations will be presented and dis- 
cussed. Finally, we will draw conclusions emphasizing the im- 
plications of our results for SN la models. 

2. The cellular burning regime and its numerical 
simulation 

Since our study aims at the flame evolution around the Gibson 
scale, which is expected to be well-separated from the scale 
of the flame width for fuel densities above 10 7 gem -3 , we will 
treat the flame in the discontinuity approximation in the follow- 
ing. That is, we describe the flame as a simple discontinuity 
between fuel and ashes and ignore its internal structure. Note 
that this picture does not account for the microphysical trans- 
port phenomena and therefore the laminar burning speed of the 
flame (i.e. the velocity of a planar deflagration flame) is not 
intrinsically given in this model. 

The Gibson scale is determined by this laminar burning 
speed si, by the velocity fluctuations at the integral scale 
of turbulence, and by the scaling law of the eddy cascade. 
iTimmes & Wooslevl(ll992l) determined the laminar burning ve- 
locity by means of one-dimensional simulations fully resolv- 
ing the internal structure of the flame, and we apply the fitting 
formula given there to obtain the value for a particular com- 
posi tion and density of the fuel. As has been mentioned by 
Roo ke et aD d2003ah . this formula is afflicted with some un- 
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certainties, which introduce a considerable ambiguity to the 
exact value of the Gibson scale. Following the approach of 
Rop ke etaD J2003al) . we will arrange our simulations at length 
scales around 10 4 cm as a crude estimate of the Gibson length, 
since our goal is to qualitatively explore the possible effects 
here rather than to perform highly precise measurements. 

Below the Gibson scale, flame propagation is determined 
by the competition between the Landau-Darrieus (LD) insta- 
bility and its nonlinear stabilization. The LD instability is of 
pure hydrodynamical origin and results from the density con- 
trast across the flame front in combination with mass flux con- 
servation. This leads to a refraction of the stream lines at the 
flame causing a flow field that enhances initial perturbations of 
the planar flame shape. The fact that the LD instability acts on 
thermonuclear flames under conditions of SN la explosions has 
been proven by means of a full hydrodynamical simulation by 
Nie mever & H illebrand d 111 99511 . The refin ed numerical model 
that was described by Ropke et al. ( 2003b) allowed us to inves- 
tigate this instability in detail. We simulated the propagation 
of a sinusoidally perturbed flame front into quiescent fuel. The 
initial linear regime of flame evolution was shown to be in- 
deed determined by the LD instability. The flow field account- 
ing for this effect was apparent in our simulations ( R opke etall 
l2003al) . Moreover, the growth rate of the amplitude of the per- 
turbation was found to be consistent with the dispersion rela- 
tion that lLand'aul ll 19441) derived by means of a linear stability 
analysis. 

The nonlinear stabilization of the flame is a geometrical ef- 
fect lZel'dovichl ill 9661) explained the mechanism by following 
the propagation of a perturbed flame by means of Huy gens' 
princip le of geometrical optics (see also Fig. 1 in iRopke et alJ 
2003b). Once the perturbation has grown to a critical size, for- 
mer recesses of the flame front develop into cusps between 
bulges of the front. It can easily be shown that the propaga- 
tion velocity v cusp of such a cusp exceeds the laminar burning 
velocity Si of the other parts of the flame by 



K cusp 



cos 6 



(1) 



where 9 denotes the inclination angle of the bulges adjacent to 
the cusp. This effect balances further growth of the perturba- 
tion due to the LD instability and causes a cellular shape of 
the flame front. Therefore the regime of flame propagation at 
scales where it is dominated by the interplay between the LD 
instability and its nonlinear stabilization will be termed cellular 
burning regime in the following. 

The nonlinear stage of flame evolution could not be 
reached by early atte mpts of hyd rodynamic al sim ulations of 
flame propagation (Nie mever & H illebrandd 1 19951) . However, 
applying a semi-ana l ytical description of flame evolution, 
Blinniko v & Sasorovl dl996l) explored effects of the cellular 
burning regi me in the conte xt of SN la explosions and also 
iBvchkov & Libermanl \ 19951) discussed the cellular flame sta- 
bilization. With he lp of th e numerical implementation that was 
described by Ropke et al. (2003b) on the basis of the work by 
iReinecke et alJ i 19991) . it was possible to prove by means of 
a hydrodynamical simulation, that cellular stabilization holds 
for thermonuclear flames in SN la explosions ( Ropke et all 



l2003blal) . In connection with flame interaction with vortical 
flows, which will be studied in the following, the most impor- 
tant results of the investigation of the cellular stabilization of 
the flame propagating into quiescent fuel were the following: 

- The flame stabilizes in a cellular shape for fuel densities 
down to 10 7 g crrT 3 . However, at the lowest fuel densities, 
the stabilization can only be revealed with high numerical 
resolution. This may be interpreted as a hint to increased 
sensitivity of the cellular flame stabilization to (numerical) 
noise at low fuel densities. 

- The evolution of the flame shape tends to a steady state con- 
sisting of only one cell filling the entire computational do- 
main. This is consistent with results from semi-analytical 
and analytical studies b a sed on the Sivashinsky-eq uation 
JGutman & Sivashinskvt Il990t iThual et all Il985l) . The 
single-cell solution establishes independently of the wave- 
length of the initial flame perturbation via merging of small 
cells. 

- Applying a simulation setup of an on average planar 
flame configuration in a finite computational domain leads 
to a dependence of the alignment of the finally emerg- 
ing domain-filling single-cell structure on the boundary 
conditions imposed transverse to the direction of flame 
propagation. In contrast to periodic boundaries, the crest 
will align in the center for reflecting boundary con- 
ditio ns. This is consistent w ith semi-analytical results 
JGutman & Sivashiiiskylll990l) . 

- The emerging cellular shape of the flame leads to an in- 
creased flame surface and equivalently an increased effec- 
tive flame propagation speed as compared to a planar flame 
front. 

- Depending on the width of the computational domain, the 
fundamental single-cell structure may be superimposed by 
a short-wavelength cellular pattern. However, these small 
cells are advected towards the cusp where they disappear. 
This mechanis m is consi stent with theoretical predictions 
hvlZerdovichet^l(ll98(i 

These results will be helpful in the interpretation of the simu- 
lations of flame propagation in vortical flows. The possibility 
of a weaker cellular stabilization at low fuel densities is a fur- 
ther motivation to test the interaction of the flame with physical 
noise, represented by a vortical flow field. 

We will only briefly outline t he num eric al implementa- 
tion he re and refer to Ropke etaD J2003bl) and lReinecke et all 
( 1999) for detailed discussions. The description of the flame 
in the discontinuity approximation allows us to separate the 
flame propagation from the hydrodynamics in a first step. The 
hydrodynamics part of the numerical implementation thus re- 
duces to the solution of the reactive Euler equations, which 
is performed applying a finite-volume approach. In particular 
w e employ the piecewise pa rabolic method (PPM) suggested 
bv lColella & Woo dward! (Il98 4) in the Prometheus implemen- 
1 etall J 



tation bv iFrvxell etalT ( il989T) . Flam e propagation is mod- 
eled using the level - set method (e.g. tosher &Sethiadll 988: 
[ Sussman et~ai]ll994 ISmilianovski et alJll997l IReinecke etall 



1999). For the correct reproduction of hydrodynamical effects, 
such as the LD instability, a precise coupling between flame 
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Fig. 1. Vortical flow field as applied in the simulations (here for 
the case p = 5 x 10 8 gcnT 3 ). The frame of reference is co- 
moving with the vortices. The snapshot is taken after 1.5 flame 
crossing times over the entire computational domain. 



and flow turned out to be essential (Rook eet alll2003bl) . This 
is provided by applying the i n-cell reconstruction/fl ux-splitting 
technique introduced bv lSmilianovski et alJ Jl997l) . 

3. Flame interaction with a vortical flow 

In the following we will present the results from simulations 
of the interaction of cellular flames with vortical flows. After 
discussing the simulation setup we will make some general re- 
marks on what can be expected in the simulations. This will 
be followed by a survey of flame evolution depending on the 
strength of the vortices for an exemplary case at a fixed fuel 
density. Finally, we will capture the simulation results for vari- 
ous fuel densities in a quantitative way. 

3.1. Simulation setup 

The basic setup of our simulations was similar t o that used to 
describe flame propagation into quiescent fuel jRopke et all 
2003b a). The flame was initialized in a computational domain 
with periodic boundary conditions transverse to the direction 
of flame propagation. In order to capture the full flame front 
in cases of strong deformation, we changed from a quadratic 
domain to a rectangle of 300 x 200 computational cells. The 
cell width was set to Ax — Ay — 50 cm. Several experiments 
with initially planar flames yielded drastic responses of the 
flame shape when the vortices encountered it. This can be at- 
tributed to the high sensitivity of planar flames to perturbations. 
Consequently, in order to make the results of different simula- 
tions comparable and to enable a quantification of the flame 
evolution, the flame was perturbed initially in a sinusoidal way 
with eight periods fitting into the domain. This initial perturba- 
tion grows due to the LD instability and stabilizes in a cellular 
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Fig. 2. Profile of the y-component of the vortical flow field par- 
allel to the x-axis (a) p = 5 x 10 8 g cm 4 ; (b) p = 1 x 10 s g cm" 3 . 
The snapshots are taken after 1.5 flame crossing times over the 
entire computational domain. 



pattern, before the injected vortices reach the flame front. In 
this way it was possible to study the interaction of a stabilized 
flame with a vortical flow field preventing the incoming flow 
from unpredictably deforming the flame shape at the first en- 
counter. 

The numerical investigation of flame interaction with tur- 
bulence is generally an intricate issue. In principle, the straight 
forward way to go would be to produce an isotropic turbulence 
field by external forcing and to set up a flame in this field. This 
approach is, however, too expensive for the purpose of the pa- 
rameter study we are aiming at. For this reason we simply mod- 
ified the inflow boundary condition in order to inject a vortical 
flow field instead of quiescent fuel. The flow field is not stirred 
and turbulence is not actively produced inside the computa- 
tional domain by external forcing. As a consequence turbulence 
injected here partly decays before reaching the flame front — an 
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Label 


p a [ gem 3 ] 


p 6 [gcm 3 ] 


P 


At 


Jilcms" 1 ] (TW) 


e,.„[ergg '] 


e,,ft[ergg '] 


1.25 x 10 7 


1.232 x 10 7 


3.830 x 10 s 


3.22 


0.526 


2.86 x 10 5 


2.005 x 10 17 


8.99 x 10 17 


2.5 x 10 7 


2.484 x 10 7 


8.69 x 10 6 


2.86 


0.482 


5.01 x 10 5 


2.56 x 10 17 


9.50 x 10 17 


5 x 10 7 


4.988 x 10 7 


2.071 x 10 7 


2.41 


0.413 


8.74 x 10 5 


3.584 x 10 17 


1.051 x 10 18 


7.5 x 10 7 


7.50 x 10 7 


3.345 x 10 7 


2.24 


0.383 


1.21 x 10 6 


4.29 x 10 17 


1.13 x 10 18 



effect which we have to take into account in the measurement 
of characteristic qua ntities. 

As SUgg ested by Helen brook & Lawldl999l) . we applied an 
oscillating inflow boundary condition on the right hand side of 
the computational domain, which generates a vortical velocity 
field approaching the flame. The velocity at the boundary now 
reads: 

Vx = si {-1 + v' sin 2fc?ry cos 2fc7r(x - ts/)} (2) 

v y — si V cos 2knys\n2kn(x — ts/). (3) 

The parameter V characterizes the strength of the imprinted 
velocity fluctuations and k denotes the wavenumber of the os- 
cillation. This produces w hat is termed "square vortices" by 
Helenbrook & La^ i 19991) . A similar flow field is applied by 
Vladimirova et al. (2003) who name it "cellular flow" and also 
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flame propagation in it. The particular flow has the advantage of 
being very simple — ideally containing only one Fourier mode. 
This helps to make the effects visible more clearly and sim- 
plifies the interpretation. Figure ^ illustrates the vortical flow 
applied in our simulations. It shows an example at a density of 
p — 5 x 10 8 g crrT 3 . In Fig. [2^ the y-component of the veloc- 
ity is plotted against the x-coordinate for the same simulation. 
Obviously, the velocity fluctuation is damped while propagat- 
ing from the inflow boundary into the computational domain. 
The damping is even stronger for lower densities. Figure [2t> 
shows a plot for p = 1 x 10 8 g cm -3 . 

In order to be able to follow the long-term flame evolution, 
the simulations are carried out in a frame of reference comov- 
ing with the flame. For a planar flame this could be achieved 
by applying an inflow condition on the right hand side of the 
domain, where fuel enters at the laminar burning velocity, and 
an outflow condition on the opposite side of the domain. Since 
flame wrinkling and subsequent surface enhancement due to 
the LD instability and interaction with turbulent eddies will 
accelerate the flame, additional measures have to be taken in 
order to keep the flame centered in the do main. In the simula- 
tions presented by Ronk e et alJ J2003blal) the grid was simply 
shifted according to the flame motion relative to s/. However, 
this approach is not consistent with the oscillating boundary 
condition. Therefore we modified eq. to 

v x — si I — h v' sin 2kny cos 2kn(x - ts/) 1 , (4) 



^ A planar J 

where A denotes the actual flame surface (measured in each 
time step) and A p i anar stands for the surface of the correspond- 
ing planar flame, which is equivalent to the width of the domain 



in our case. This ensures that the velocity with which the fuel 
enters the domain adapts to the current effective flame propaga- 
tion speed (cf. eq. (jfji in Sect. 13. 4> and the position of the flame 
relative to the grid is kept approximately fixed. 

We performed a number of simulations with different fuel 
densities. For convenience of the reader, a compilation of the 
setup values is given in Tabled 

3.2. What can be expected? 

Before we present the results of our numerical study we discuss 
what we actually are looking for in the simulations. That is, we 
somehow have to define how to discriminate between stable 
flame propagation and a breakdown of the stabilization. Why 
is this a question of definition? From the conjecture of active 
turbulent combustion it may appear that a breakdown of the sta- 
bilization must be obvious from a rapid and unlimited growth 
of the flame propagation velocity. However, this cannot be ex- 
pected to occur in a numerical simulation where the growth of 
the flame surface due to ATC would be limited by discretization 
and resolution. Hence, the maximum effect that can be antici- 
pated is a rapid increase in flame surface and burning velocity 
and a saturation at some higher value. A second possibility is 
that with increasing strength the vortical flow starts to dominate 
the flame evolution. 

On the other hand, how can we guarantee flame stabiliza- 
tion? Owing to the high computational costs it will not be possi- 
ble to follow the flame propagation for an arbitrarily long time. 
However, in case of stabilization we can actually make use of 
the knowledge of the flame pattern that is to be expected in case 
of stability in our specific simulation setup. From the prece ding 
studies o f flame propagation into quiescent fuel ( Ronke et all 
2003b a) it is known that the flame finally stabilizes in a single 
domain-filling cusp-like structure for our setup. This structure 
may be superimposed by a smaller-scale cellular pattern in suf- 
ficiently resolved simulations. A similar result can be antici- 
pated for interaction with weak imprinted vortices. 

Thus we may expect two extremal behaviors of the flame. 
In one case the flame stabilizes and the incoming perturbation 
fails to break up the cellular pattern. The cells of the stabilized 
small wavelength pattern caused by the initial perturbation of 
the flame front will then merge forming the single domain- 
filling cusp which then propagates stably. Alternatively, if the 
intensity of the incoming vortices is high enough to destroy the 
stabilization, no single domain-filling cusp will finally emerge. 
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Fig. 4. Flame evolution for a fuel density of 5 x 10 7 gcirT 3 and a velocity fluctuation of v'/si = 2.5 at the right boundary. 



Snapshots taken at time steps (a) t = 4.0 x 1CT J s, (b) t = 8.0 x 10~ j s, (c) t = 3.2 X 10~ 2 s, and (d) f = 4.8 X 10 _/ s. The vorticity 
( (cf. eq. Q) is color-coded and the flame position is indicated by solid white curves. 



In this second case the flame should rather show a transient 
pattern. 

3.3. General features 

As an example, we consider the case of a fuel density amount- 
ing to 5 x 10 7 g cirT 3 (cf. Table^. Figure[3]gives an overview 
over the flame evolution in interaction with vortical fuel flows 
of different strengths. The values v' /si = 0.7, v'/si = 2.0, and 
v'/si = 2.5 indicate the amplitudes of velocity fluctuations im- 
posed at the right boundary of the domain. Note, however, that 
they do not necessarily represent the values experienced by the 
flame front for reasons given in section l3~Tl 

Figure [3^ provides the comparison with flame propaga- 
tion into quiescent fuel for the chosen setup. In this simula- 
tion v' was set to zero. It resembl e s the fe atures of the sim- 
ulations presented by Ron ke etahl J2Q03hlah . Due to the LD 
instability, the initial perturbation grows and the flame shape 
evolves into a cellular pattern in the nonlinear regime. The sub- 
sequent contours show the "merging" of the short-wavelength 
cells imprinted by the initial condition resulting in the for- 
mation of larger cells. This behavior was observed in a sim- 



ulati on of the long -term flame evolution in quiescent fuel 
(Ro pke et alll2003al) . There a cusp-like structure that was cen- 
tered in the domain finally emerged as the steady state struc- 
ture. The tendency to form a domain-filling cell is also apparent 
in Fig. [3^. However, following the evolution until the steady- 
state is reached is very expensive and shall not be tried here. 
This is certainly a drawback in the current study, but it was 
chosen as a compromise in order to be able to explore a larger 
parameter space with given computational resources. 

Contrary to the flame evolution in Fig. [3^, the case of 
v'/si = 2.5 (Fig.|5Jl) clearly shows the disruption of the ini- 
tial cellular pattern by interaction of the flame with vortices. 
Also, the formation of a single domain-filling cell is suppressed 
by the interaction. The overall flame shape evolution in this 
case can be interpreted as an adaptation to the imprinted vor- 
tical flow structure. This is emphasized by Fig|4] It provides a 
more detailed visualization of the interaction between the flame 
shape and the vortical flow, which is characterized by the color- 
coded logarithm of the absolute value of the vorticity ( of the 
flow field, 



dv y dv x 
dx dy ' 



(5) 
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Fig. 5. Flame evolution for a fuel density of 5 x 10 7 gcm 3 and a velocity fluctuation of v'/si = 0.7 at the right boundary. 



Snapshots taken at time steps (a) t = 8.0 x 10" J s, (b) t = 2.4 x 10~ z s, (c) t = 3.2 X 10~ 2 s, and (d) t = 4.8 X 10 _/ s. The vorticity 
( (cf. eq. Q) is color-coded and the flame position is indicated by solid white curves. 



pointing out the vortices imprinted on the fuel flow. 

The transition in the flame evolution between the extreme 
cases depicted in Figs. [2t and 01 does, however, not proceed 
abruptly, but rather in a smooth transition. FiguresOJ>,c present 
two examples of intermediate behaviors at v'/si - 0.7 and 
v'/si = 2.0, the first of which is again illustrated by snapshots in 
Fig-0 In this case, still a tendency to form a long-wavelength 
cellular structure is visible. Here the mechan ism of the advec- 
tion of small perturbations t oward the cusp dZel'dovich et all 
Il980t iRopke et all l2003blal) dominates over the effects of the 
vortical flow on the flame shape. 

Simulations similar to the aforementioned were also per- 
formed in some cases with higher resolution of the numer- 
ical grid in order to test whether or not insufficient reso- 
lution prevents small scale flame structures from develop- 
ing. However, no significant difference to the flame evolu- 
tion depicted in Fig . Et-d could be noticed. This agrees with 
iHelenbrook & Law! dl999l) . who found in their simulations of 
chemical flames that the wavelength of perturbations that de- 
velop in the interaction with a vortical flow is determined by 
the scale of the vortices and not by a critical wavelength (in- 



troduced by a curvature-dependent burning speed of the flame 
front according to Markstein 1951). 



3.4. A Parameter study 

In order to quantify the flame evolution as qualitatively de- 
scribed in the preceding section, we conducted a parameter 
study. It aimed at the determination of flame behavior as a func- 
tion of the strength of the vortices in the incoming flow and in 
dependence on the fuel density. 

As discussed by .'Rop ke et"aD d2003al) . the simulation of 
flame stabilization at around 1.0 x 10 7 gcm~ 3 requires unaf- 
fordable high numerical resolutions and furthermore our thin 
flame approximation is certainly not a realistic description of 
thermonuclear burning in SNe la at these low densities. On 
the other hand, it can be observed, that at higher fuel densities 
(around 10 8 g cm 4 ) the sensitivity of the initial flame to small- 
wavelength perturbations increases. This causes difficulties in 
producing a stable flame configuration prior to the incoming 
vortices encountering the flame front. However, in connection 
with a possible DDT we are mainly interested in the late stages 
of the SN la explosions (see Sect. Therefore we restricted 
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Fig. 3. (a) Flame propagation into quiescent fuel at a density 
of 5 x 10 7 gcirT 3 and interaction with a vortical fuel field for 
velocity fluctuations at the right boundary of (b) v'/s/ = 0.7, 
(c) v'/si = 2.0, and (d) v'/si = 2.5. Each contour represents a 
time step of 3.2 x 10~ 3 s. 

the study to four fuel densities: p u = 1.25 x 10 7 gctrT 3 , p u = 
2.5 X 10 7 g cnr 3 , p u = 5 x 10 7 g cnr 3 , and p u = 7.5 x 10 7 g cm 3 
(cf. TableQ). 

Plots similar to Fig. Q] correspon ding to simu lations with 
different fuel densities are given by|R6pke (2003). The flame 
evolution in these simulations can be summarized as follows: 
At lower fuel densities the effect of the incoming vortices 
on the flame structure is generally more drastic. For p„ = 
2.5 x 10 7 gem -3 the tendency of the formation of the domain- 
filling structure is still visible in the case of propagation into 
quiescent fuel. This is not the case for p„ = 1.25 x 10 7 gem -3 



for reasons discussed by Ropke etafl d2003al) . However, with 
increasing strengths of the imprinted velocity fluctuations, the 
flame still gradually adapts to the flow for both fuel densities. 
At p u — 7.5 x 10 7 g crrT 3 the flame evolution is similar to that 
at p u — 5 x 10 7 g cirT 3 . Only in case of propagation into quies- 
cent fuel an initial destabilizatio n with resp ect to small wave- 
length perturbations (cf. IRopke et alJl2003al) alters the evolu- 
tion for a transition period, before the flame finally stabilizes in 
a domain-filling single-cell structure. 

In a quantitative evaluation of our parameter study regard- 
ing the flame evolution for different fuel densities we will now 
address (i) the dependency of the effective flame propagation 
velocity on the strength of the imprinted velocity fluctuations 
and (ii) the amplification of the velocity fluctuation by the 
flame front. 

For the exemplary case of p„ — 5 x 10 7 gcirT 3 we discuss 
the measurement of the necessary quantities in our simulations. 
The effective propagation velocity v e ff of the flame is deter- 
mined via the flame surface: 



v eff _ A sur f 
$1 ^planar 



(6) 



Underlying to that relation is the assumption of a constant, ge- 
ometry independent burning velocity s/ of the flame. This as- 
sumption has been disc u ssed in the context of the present flame 
model by IRopke et alJ J2003bl) . The measured flame surface 
area as a function of time is plotted in Figs. [(J and|6j? for the 
simulations with v'/si - 0.7 and v'/si = 2.5, respectively. Here 
adaptation of the flame structure to the strong vortical fuel flow 
is obvious. After the imprinted velocity fluctuations reach the 
flame at t ~ 150t Av , its surface becomes dominated by the vor- 
tices and changes with the frequency of these. This is reflected 
by the high-frequency fluctuations of A in Fig. [(J). In addition 
to fluctuations on a small time scale, the interaction with the 
vortical flow field introduces considerable long term fluctua- 
tions. These are, of course, expected in the initial phase when 
the flame begins to react to the incoming vortices. But even at 
much later times the flame surface area does not reach a steady 
state. 

The standard deviation of the velocity field is obtained sep- 
arately in the fuel upstream of the front and in the ashes down- 
stream of it in order to determine the amplification of the ve- 
locity fluctuations across the flame front. In Fig.Qthe quanti- 
ties are plotted against time for our exemplary case. As men- 
tioned above, the strength of the vortices produced at the inflow 
boundary is not a reliable measure of what the flame actually 
experiences, since the velocity fluctuations get damped quickly 
when propagating toward the flame and bending of the flame 
front may lead to different velocity fluctuations at different lo- 
cations on the flame. Moreover, the dependence of the damping 
on the fuel density hampers the comparability of simulations 
with this parameter varying. Therefore we determine the stan- 
dard deviations of the velocity fields, cr(v„) and cr(vi), within 
a certain belt around the flame front. This belt was scaled in a 
way that the time for flame crossing over its width remained 
the same for all fuel densities. 

All quantities measured in our simulations fluctuate consid- 
erably with time. Hence we averaged them over a time ranging 
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Fig. 7. (a,b) standard deviation of the velocity field ahead of the front over time, and (c,d) ratio of the standard deviations of 
the velocity fields beyond and ahead of the flame over time. The left and right columns of plots correspond to v'/si = 0.7 and 
v'/s/ = 2.5, respectively. 



from 480 ta.i to 720 ta. v (Ta.i denotes the crossing time of the 
corresponding laminar flame over one grid cell). 

In Fig. [8] the temporal mean value of flame surface area 
A sur f (normalized to the surface area of the corresponding pla- 
nar flame front A p i anar ) is plotted against the temporal mean of 
standard deviation of the velocity in the fuel region (normal- 
ized to the laminar burning velocity of the flame). The increase 
of the flame speed with stronger incoming velocity fluctuations 
is evident here. However, there is quite a large scatter in the 
data. This has two reasons. First, we certainly did not follow the 
flame evolution long enough to undoubtly reach the steady state 
of the flame shape — which should ideally be the clearly dis- 
tinct single domain-filling cusp-like structure in case of weak 
incoming vortices. To follow the flame propagation up to this 
stage would, however, be much too expensive for a parameter 
study. Second, some peculiar features of the flame shape may 
develop as a result of small perturbations (e.g. due to numeri- 
cal noise). In particular, this may cause substantial deviations 
if it affects the long-wavelength structure of the flame and is 
primarily responsible for the scatter of the results in case of 



strong incoming velocity fluctuations. This scatter could very 
lik ely be cured by taki ng the mean over longer time intervals. 

lDamkohlerl ll940) proposed a linear dependence of the ef- 
fective flame propagation velocity on the turbulence intensity 
for turb ulent c ombustion in what is today called the fiamelet 
regime JPetersl [l986). Our results are consistent with a linear 
growth law. Corresponding fits to the data are included in Fig. [3] 
and values for the parameter a according to the fitting formula 



<A SU rf> , „ <0-(v„)> 

= 1.0 + a 

Apianar $1 



(7) 



are given in Table |2] It is evident from Fig.[8]that interaction of 
the flame in the cellular regime with turbulent velocity fluctu- 
ation can lead to a substantial acceleration of the flame propa- 
gation velocity. Unfortunately, it is not possible to infer a clear 
trend of the slope depending on the fuel density from this pa- 
rameter study. 

In order to estimate the amplification of the velocity fluc- 
tuations across the flame for different densities, we compared 
the strength of velocity fluctuations downstream of the flame 
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Table 2. Fit parameters according to fitting formulas Q and © with corresponding asymptotic standard errors. 



p„ [gem 3 ] a b c 

1.25 x 10 7 1.25712 ±0.0215 0.177654 ± 0.05236 1.28103 ± 0.0.04581 

2.5 x 10 7 1.29816 ±0.0797 0.180935 ± 0.01718 0.779895 ± 0.01982 

5 x 10 7 1.05533 ±0.04931 0.141428 ± 0.01171 0.558411 ± 0.01883 

7.5 x 10 7 1.07245 ±0.05131 0.150079 ± 0.2911 0.466975 ± 0.03134 



c 
I 




200 400 600 800 




200 400 600 800 



Fig. 6. Flame surface area over time (normalized to the crossing 
time of the laminar flame over one grid cell tax) for (a) v'/si = 
0.7 and (b) v'/si = 2.5. 

with the strength of the imprinted vortices in the fuel. As al- 
ready shoyvTi_^orthe_case of flame propagation into quiescent 
fuel dRopke et all l2003al) . the flame produces vorticity in the 
ashes. Here, the question is addressed whether vorticity present 
in the fuel will be amplified in the flame. 




0.0 0.5 1.0 1.5 

(cr(v u ))/si 



Fig. 8. Dependence of the flame surface area on the strength of 
the imprinted velocity fluctuations. 

The result of the study for a variety of fuel densities is 
plotted in Figure |9] In this plot, trends show up much more 
clearly than in Fig. [8] However, some scatter is still present in 
the data. As can be seen from snapshots of the evolution of the 
flame front (an example is the bottom left snapshot in Fig.|4j, 
the merging of cusps can eventually produce transient "bursts" 
in the vorticity downstream of the flame. This is particularly 
prominent in case of changes in the long-wavelength flame 
structure, but not necessarily connected to it. From Fig. |7Jl, 
where this "burst" in velocity fluctuation appears as a spike in 
the profile, it can be concluded that the duration of those events 
is very short. Hence their contribution to the temporal mean is 
small (whereas the long-wavelength flame structure and thus 
the flame surface changes slowly). Consequently, the scatter in 
plot[9]is much smaller than in plot|8] 

The data plotted in Fig. [9] can be fit rather well by a lin- 
ear relation between the ratio (cr(vi)/cr(v„)) and the incoming 
turbulence intensity, i.e. 

(<r(v b )) = b + c(oiv u )). (8) 

Table|2]provides the fit parameters. In contrast to Figure[8] the 
plot in Figure [9]reveals a clear trend for the fuel density. The 
slope c increases with lower p u . This is what would be expected 
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Fig. 9. Relation between velocity fluctuations upstream and 
downstream of the flame front. 
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Fig. 10. Ratio of the velocities upstream and downstream of the 
flame front as a function of the density contrast fi over it. 



taking into account the increased production of specific volume 
across the flame with lower p„. In Fig.^|the factor c is plot- 
ted against the corresponding density ratio fi, accentuating the 
trend. Of course, a functional dependence cannot be inferred 
from this small sample of data. It is interesting to note that for 
p u ~ 2.5 x 10 7 gcrrT 3 the slope c becomes greater than unity, 
indicating a slight amplification of velocity fluctuation in this 
case. 



The above interpretation of this result requires some cau- 
tion, since the cellular stabilization of the flame has been shown 
to be weaker with lower fuel density. However, this seems to 
be a numerical rather than a physical effect, since the cusps be- 
come more stable with higher resolution. Thus, at the given res- 
olution of our parameter study, the flame evolution shows pe- 
culiar behavior below a fuel density of 1 .25 x 1 7 g cirT 3 . Here, 
cusps may loose stability distorting the flame shape, which then 
responds with accelerated propagation of the perturbed part for 
a short period, after that the flame stabilizes again. However, 
we do not observe such a flame evolution at higher fuel den- 
sities included in this parameter study. Thus, the dependency 
on the fuel density is likely to be (at least partially) of physical 
origin. 

4. Conclusions 

In the presented simulations we performed the first attempt to 
investigate the interaction of cellular flames with turbulent ve- 
locity fluctuations in the context of SN la explosions. The re- 
sults can be summarized as follows: 

1 . In case of the interaction of the flame with weak vortices 
in the fuel flow, the flame exhibits an evolution similar to 
that observed in case of flame propagation into quiescent 
fuel. In the setup chosen for our simulations the flame front 
showed a tendency to align in a single domain-filling cellu- 
lar structure. Interaction with the incoming vortices leads to 
a superposition of this fundamental structure with a cellu- 
lar pattern of shorte r wave l ength. S imilar to the simulations 
presented by Roo ke etafl d2003bl) . where a cellular super- 
position of the basic flame structure was observed in case of 
wide computational domains, this does not lead to a desta- 
bilization of the domain-filling cel l . In a ccord with theoret- 
ical predictions dZeFdovich e t al.. 198(J) the small super- 
imposed perturbations are advected toward the cusp where 
they disap pear. The flow field accounting for this effect was 
studied by Ro pke et all d2003al) . 

2. If the cellularly stabilized flame interacts with vortices of 
sufficient strength, then the initial cellular pattern will break 
up in the sense that the initial shape is lost and the flame 
adapts to the incoming vortical flow. The tendency of the 
flame to form a single-cell structure is suppressed by this 
effect. However, we do not observe a sharp transition but 
rather a gradual change in the flame behavior. 

3. Our numerical simulations suggest that the flame smoothly 
adapts to the incoming vortical flow with increasing turbu- 
lence intensity. No drastic effects were observed. Although 
we measured an amplification of the strength of the veloc- 
ity fluctuation across the flame forp„ = 2.5 x 10 7 gcm -3 , 
we could not observe a deviation from a linear scaling be- 
tween the flame surface area (and thus the effective flame 
propagation speed) and the strength of the imprinted veloc- 
ity fluctuations. 

We draw the conclusion that the results of our study are 
consistent with flame propagation in the flamelet regime of 
turbulent combustion, at least at fuel densities above p u » 
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2.5 x 10 7 gcrrr 3 . Thus, our study corroborates the assump- 
tion of flame stability at unresolved scales in large-scale SN 
la explosion models. This is an important contribution that our 
small-scale model makes to the credibility of those simulations. 
The slight amplification of the turbulent velocity fluctuations 
across the flame at p„ = 2.5 x 10 7 gcirT 3 (and possibly even 
more pronounced below that value) may have no significant 
impact on the current SN la models, since the flame width at 
those densities will become n on-negligible (cf. the data given 
bv lTimmes & Wooslevl 19921 . Then the burning is expected to 
enter the distributed regime of turbulent combustion, whose 
simulation requires completely different numerical methods. 
In the context of SN la explosions, an approa ch to study this 
regime was discussed bv lLisewski et alJ (2000). 

Although the increase in effective flame propagation result- 
ing from the cellular regime is negligible compared to the in- 
crease in flame speed at larger scales where the flame is af- 
fected by the turbulent cascade, it may be significant in the 
early sta ges of the explosion . In the current large-scale mod- 
els (e.g. Rei necke et al.ll2002l) . the flame is assumed to propa- 
gate with its laminar flame velocity before the rising Rayleigh- 
Taylor bubbles establish the turbulent cascade. This may cause 
a proble m with the nucl eosynthetic yields from the SN la ex- 
plosion ( Rein ertell2001 ). Since the laminar flame speed is very 
low, the WD expands slowly in the beginning of the explo- 
sion and thus the combustion products remain at rather high 
densities for a considerable time. This effect causes a neu- 
tronization of the material by electron capt ure and an overpro- 
duction of neutron-rich heavy nuclei ( Nom oto & Kon do , 1 99 1 ; 
Brac hwitz et all l2000l) . eventually even leading to a collapse 
of the star. Taking into account the velocity increase result- 
ing from the cellular regime, this problem could be extenu- 
ated. As has been shown in this study, this may be in partic- 
ular the case if strong velocity fluctuations left over from the 
pre-ignition convection interact with the flame. However, the 
strength of turbulenc e resulting from pre-ign ition effects is not 
well-determined yet. iHoflich & Steinl J2002t) claim values as 
high as ~10 7 cm s _1 . The impact of the cellular burning regime 
on the nucleosynthesis yields of the SN la explosion models 
will be investigated in a forthcoming study. 

Even though we can probably not completely rule out the 
possibility of active turbulent combustion, we found no con- 
vincing hint for such an effect in our numerical investigations. 
Our simulations indicate that effects resulting from the cellu- 
lar regime of flame propagation are unlikely to trigger a pre- 
sumed deflagration-to-detonation transition. Hence, the search 
for active turbulent combustion and deflagration-to-detonation 
transition should focus on the distributed burning regime. 
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